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We investigate the possibility of an Anderson type transition in the quantum kicked rotor with 
a smooth potential due to dynamical localization of the wavefunctions. Our results show the typ- 
ical characteristics of a critical behavior i.e multifractal eigenfunctions and a scale-invariant level- 
statistics at a critical kicking strength which classically corresponds to a mixed regime. This indicates 
the existence of a localization to derealization transition in the quantum kicked rotor. 

Our study also reveals the possibility of other type of transitions in the quantum kicked rotor, with 
a kicking strength well within strongly chaotic regime. These transitions, driven by the breaking 
of exact symmetries e.g. time-reversal and parity, are similar to weak-localization transitions in 
disordered metals. 
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INTRODUCTION 



The analogy of the statistical fluctuations of dynamical systems and disordered systems is well-known in the 
delocalized wave-regime (corresponding to the metallic limit in disordered systems and the classically chaotic limit 
in dynamical systems) and has been explained using random matrix theory as a tool [3, [2, [H, IH • A similar analogy 
exists for fully localized regimes of the wavefunctions too (i.e. between the insulator limit of disordered systems 
and the integrable limit of dynamical systems) [l|, 0, 9- It is therefore natural to probe the presence/ absence 
of the analogy in partially localized/ critical regimes of these systems. Our analysis shows that, similar to the 
d > 2 Anderson Hamiltonian (d as dimension), the d = 1 quantum kicked rotor (QKR) undergoes a localization- 
delocalization transition in the classically mixed regime. We also find quantum phase transitions in its chaotic regime 
due to breaking of the symmetries e.g. time-reversal and parity in the quantum system. Similar to disordered systems, 
the symmetry-breaking transitions in the QKR occur due to weak-localization effects. Similar phase transitions due 
to symmetry-breakings have been seen in a few other complex systems too e.g. the ensembles of distinguishable spins 

The connection of the kicked rotor to the d = 1 Anderson Hamiltonian has been known for several decades 
@> B H HI- A recent work Q further explores the connection and shows that, for the non-analytic potentials in 
the QKR, the eigenstates show multifractality or power-law localization [l(| [Hi G2> EH , a behavior similar to the 
eigenstates of a d > 2 dimensional Anderson Hamiltonian at its critical point. Our study shows existence of 
the multifractal eigenstates in the QKR with smooth potentials too e.g. V(q) = Kcos(q) at specific parametric 
conditions. Furthermore, similar to a critical Anderson system, the multifractality in the QKR is accompanied by a 
critical level statistics (size- independent and different from the two ends of the transition), a necessary criteria for the 
critical behavior [To| . This indicates a much deeper connectivity of the kicked rotor to the Anderson Hamiltonian, 
not affected just by the nature of the potential or the dimension of system. As discussed here, the connection seems 
to be mainly governed by the "degree of complexity" (measured by the complexity parameter discussed later) and 
may exist among a wider range of dynamical and disordered systems. 

The present study is motivated by a recent analytical work [lH [l6| leading to a common mathematical formulation 
for the statistical fluctuations of a wide range of complex systems. The work, based on the ensemble-averaging, shows 
that the fluctuations are governed by a single parameter A besides global-constraints on the system [l5|, [l|| . A referred 
as the complexity parameter, turns out to be a function of the average accuracy of the matrix elements, measured in 
units of the mean-level spacing. The fluctuations in two different systems, subjected to similar global-constraints, are 
analogous if their complexity parameters are equal irrespective of other system-details. 

The A-formulation was recently used by us to find the Gaussian Brownian ensemble (GBE) [TtJ analog and the 
power-law random banded matrix ensemble (PRBME) [l8| analog of the Anderson system (for arbitrary d) 
However it can not directly be applied to find the QKR analog; this is, in principle, due to inapplicability of the 
ensemble averaging to dynamical systems. Fortunately it is possible to derive A for dynamical systems by a semi- 
classical route, using the phase-space averages [20]. The "semi-classical" A was used by us [l^,[2l| to map the statistics 
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of the time-evolution operator U of the QKR to the circular Brownian ensembles (CBE) [l7|. The lack of a suitable 
criteria for the critical statistical behavior prevented us earlier from a critical QKR-analysis. Our present work pursues 
the analysis by first analytically identifying the critical QKR-behavior using the "semi-classical" A = Afc r ; the limit 
N — > oo A — > A*, with A* independent of the size N, gives the critical points of transition. This is followed by a 
numerical analysis of an ensemble of the QKR which confirms critical behavior at the semi-classically predicted values. 
This in turn suggests a paradoxical validity of the ensemble averaging in dynamical systems at least of QKR type; (also 
indicated by the CBE-QKR mapping). A subsequent comparison of the "semi-classical" Afc r to the "ensemble-based" 
A ae of the Anderson Hamiltonian (arbitrary d) gives us its QKR analog; the analogy is numerically confirmed too. 

The statistical behavior in the bulk of the spectrum of a standard Gaussian ensemble is known to be analogous 
to a standard circular ensemble (for large matrix sizes) [13, [HJ . Our work extends this analogy to their Brownian 
ensemble counterparts too, that is, between the GBE and the CBE (by the mapping GBE — > AE — > QKR — > CBE). 
This in turn indicates a connection among a wide-range of physical systems which are known to be well-modeled by 
the Anderson Hamiltonian, the kicked rotor and the Brownian ensembles of Gaussian and circular type. 

The paper is organized as follows: The section II briefly reviews the basic features of the kicked rotor and the 
Anderson Hamiltonian required for our analysis; it also discusses the parametric conditions in the kicked rotor which 
can support critical points. The section III deals with the numerical confirmation of the critical level-statistics and 
the multifractality of the eigcnfunctions at critical parametric conditions in the QKR and their comparison with a 
d = 3 dimensional Anderson Hamiltonian. We conclude in section IV with a brief discussion of our main results and 
open questions. 



II. KICKED ROTOR AND ANDERSON HAMILTONIAN 



The kicked rotor and the Anderson Hamiltonian have been subjects of intense study in past and many of their 
details can be found in several references [l], 0, 0, H, [H| • However, for self consistency of the paper, we present here a 
few details required for later discussion. 

A. Kicked Rotor and Complexity Parameter 

The kicked rotor can be described as a pendulum subjected to periodic kicks (of time-period T) with Hamiltonian 
H given as 

H= {jp + ^ 2 +Kcos(9 + e ) £ S(t-nT). (1) 



n— — oo 



Here K is the stochasticity parameter, 7 and 9q are the time-reversal and parity symmetry breaking parameters in 
the quantum Hamiltonian while acting in a finite Hilbert space. 

Integration of the equations of motion = — p — ^ between subsequent kicks e.g. n and n + 1 gives the 
classical map, 

Pn+i = p n + K sin(0„ + O ) (mod 2ir), 

On+i = n + Pn+i (mod 2tt) (2) 

The map is area-preserving and invariant under the discrete translation 9^9 + 2n,p — > p + 2n. It also preserves the 
time-reversal symmetry p — > 2n — p,9 — > 9,t — > —t and the parity p — > 2it — p, 9 — > 2n — 9 for all values of 7 and 9q. 
Thus the classical dynamics depends only on K, changing from integrable (K — 0) to near integrable (0 < K < 4.5) 
to large scale chaos (K > 4.5). 

The quantum dynamics can be described by a discrete time-evolution operator U = G 1 / 2 .B.G 1 / 2 0,0] where 

B = exp [-ifccos(0 + O )] (3) 



G = exp [-iT(p + -/) 2 /2h] = exp 



(4) 
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Here k = K/h, r = Th, 8 and p — Hi-gg are the position and momentum operators respectively; p has discrete 
eigenvalues, p\m >— mh\m > (to = 1 — > N) due to periodicity of 8, (9^9 + 2ir). The choice of a rational value 
for t/Att = M/N results in a periodicity also for momentum operator p' = p + 4%M/T, (I' = I + N) and therefore in 
discrete eigenvalues for 8, (8\l >= (2Trl/N)\n >). The quantum dynamics can then be confined to a two dimensional 
torus (with a Hilbert space of finite size N = 27rmo/r with too an integer). The classical analog of this model 
corresponds to the standard mapping on a torus of size 2-nmnjT in the momentum p; thus the classical limit is 
t — > 0,k — ► oo, N — > oo with K — constant and Nt = constant @, 0j- 

For the dynamical- localization analysis, it is useful to express the matrix U in the momentum-basis 0, 0jH3] : 



1 Nl 
U mn = — exp [iT(m-7/n) 2 /4 + iT(n-7/?i) 2 /4] ^ exp [—i k cos (2irl/N + 8 )] exp [— 27ti/(to — n)/A/] (5) 

l=-Ni 

where n,m = —N±, ...,Ni with Ni — (N — l)/2 if N is odd and Aq = AT/ 2 if A" is even. It is clear that the properties 
of H, eq.JT]), are recovered in the infinite matrix size limit. 

The quantum dynamics under exact symmetry conditions (7 = 0, 80 — 0) can significantly be affected by relative 
values of k, r, N. It was first conjectured 0] and later on verified [2(| that the statistical properties are governed 
by the ratio of the localization length £ to the total number of states iV or, equivalently, k 2 /N (as £ = D/2t 2 with 
D ps K 2 /2 as the diffusion constant). However, to best of our knowledge, the critical behavior at k 2 ~ N was not 
probed before. The other parameters playing a crucial role in the quantum dynamics are 7 and 8q, the measures of 
time-reversal and parity-symmetry breaking respectively (with < 7 < h and —n/N < 8q < tt/N). Note the change 
of p — > p + 7 or 9 — > 9 + 6q is a canonical transformation, thus leaving the classical Hamiltonian unaffected. The 
corresponding quantum dynamics, however, is affected as the quantum Hamiltonian acting in a finite Hilbert space 
may not remain invariant under a unitary transformation. 

Following eq.([5]), the multi-parametric nature of U is expected to manifest itself in the statistical behavior of its 
eigenvalues (quasi-energies) and the eigenfunctions. However, as shown in [2(| using semi-classical techniques, the 
quasi-energy statistics of U is sensitive to a single parameter A^ r and the exact symmetry-conditions. Under exact 
time-reversal symmetry (7 = 0) and partially violated parity (#0 ^ 0) (taking T = 1, equivalently, r = h, without 
loss of generality), we have [2(| 

9 2 Nk 2 N*9 2 K 2 

Note, for 9q = ir/2N, Ak r ,t is essentially the same as the one conjectured in [2| for scaling behavior of the spectral 
statistics. Similarly, for the strongly chaotic case (k 2 > N) with only parity symmetry (#0 = ir/2N) and no time- 
reversal (7 ^ 0), 

7 2 A^ 3 

Akr ' nt = 48^ (7) 

with 7 g = j/h. As eq.® indicates, A-kr,t ~~ * 00 m the strongly chaotic limit k — > 00; similarly, from eq.([7]), Afc r t — * 00 
for 7 ~ 7i; the statistics in these cases can be well-modeled [20|, |2J] by the standard random matrix ensembles of 
unitary matrices, known as the standard circular ensembles e.g. circular orthogonal ensemble (COE), circular unitary 
ensemble (CUE) etc[I3]- The cases with a slow variation of k, 7 or 9q (partial localization, partial time-reversal or 
parity-violation respectively) and finite size Af correspond to a smooth variation of Ak r ,t or Kkr.nt between and 
00 . The intermediate statistics for these cases 0, [2l| can be well-described by the circular Brownian ensembles. 

The latter are the ensemble of unitary matrices, described as U w = Uq^ 2 exp[iwV]U^* and characterized by a single 



parameter A c f, e = w 2 (\Vki\ 2 ) / D 2 (D = 2n/N) and the exact system-symmetries [I7|,l22(. The perturbation V belongs 
to a standard Gaussian ensemble of the Hermitian matrices e.g Gaussian orthogonal ensemble (GOE), Gaussian 
unitary ensemble (GUE) etc [ItJ • The circular Brownian ensemble analog of a QKR is given by the condition 

Afc r = A cfce . (8) 

As mentioned above, the studies [13,[2l[ did not explore the infinite size limit of A^r and its application for the critical 
behavior analysis; we discuss it in next section. 
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B. Critical behavior of Quantum Kicked Rotor 



For a critical point analysis of the QKR, we search for the system conditions leading to the critical eigenvalue- 
statistics and the multifractal cigcnfunctions in the infinite size limit (N — > oo). The parametric conditions for the 
critical level statistics, characterized by a non-zero, finite A in the limit N — > oo, can be obtained from eqs.([6l [7J. 
The analysis suggests the possibility of several continuous families of critical points, characterized by the complexity 
parameter and the exact symmetries; we mention here only the three main cases: 

(i) k oc >/N, 7 = 0, 9q = n/2N: Eq.© in this case leads to a size-independent Ak r ,t = X 2 /16 with x — h/y/N. For 
small- A: (in the mixed regime), this corresponds to a localization — ► delocalization transition under the time-reversal 
conditions (no parity symmetry) with a continuous family of critical points characterized by \. The bulk-statistics 
here is analogous to a circular Brownian ensemble U w (see eq.®) with w = xk/zV^N (due to a GOE type perturbation 
V with |Vfcz| 2 = (1 + Ski) of a Poisson matrix [l7|, see eq.©). The two ends of the transition in this case are the 
Poisson (A — > 0) and the COE ensemble (A — > oo). 

(ii) 7 9 oc N~ 3 /' 2 ,8 — n/2N: Eq.([7j) in this case gives A kr ^ nt = \ 2 /48ir 2 with A = ^ q N 3 / 2 . For large k (in strongly 
chaotic limit), a finite A gives the critical parameter for the transition from a time-reversible to a time- irreversible 
phase (both phases delocalized) ; it can be referred as the weak-localization critical point. The end-points of the 
transition are the COE (A — * 0) and the CUE (A — * oo), with the critical statistics given by an intermediate circular 
Brownian ensemble with w — X/(2Ny/3) (due to a GUE type perturbation V, with \Vki\ 2 = 1, of a COE matrix [l7| 

(iii) 7 = 0,6»o oc N~ 3 / 2 : Here eq.® gives A kr ,t = 4> 2 K 2 / (<6An A M 2 ) with 4> = N 3 / 2 9 describing the critical point 
family for the transition from a parity-symmetric phase to a parity fully violated phase (both phases time-reversal). 
For K in the mixed regime (K < Nh), a variation of leads to the Poisson — -> COE transition. For K in the strongly 
chaotic regime (K > Nh), the transition end-points are the 2-COE (A — ► 0) and the COE (A — > oo). The critical 
statistics for a specific K is given by the intermediate Brownian ensemble with w — K4>/(4ttMN) (due to a GOE 
type perturbation V, with \Vm\ 2 = 1, of the Poisson ensemble if K << Nh, or, the 2-COE ensemble if K » Nh). 

As eqs.([51[7]) indicate, a size-independent A^ r can be obtained by other combinations of A;, 0,7 too. This suggests 
critical behavior in the symmetry-spaces other than those mentioned above. 

The critical nature of the system for specific parametric conditions can further be confirmed by an analysis of the 
eigenfunction fluctuations. The studies of a wide range of systems (see [T3, EH, El, EH and the references therein) reveal 
the presence of strong fluctuations in the eigenfunctions near a critical point. The fluctuations can be characterized 
through the set of generalized fractal dimension D q or r q — (q — l)D q , related to the scaling of the q th moment 
of the wavefunction intensity \(j){r)\ 2 with size TV [ll|: P q = J dr|</>(r)| 2<? = N- r( - q ^ d . The multifractality of the 
eigenfunctions can also be analyzed through the spectrum of singularity strengths f(a) [l(|, related to r(q) by a 
Legendre transformation: f(a(q)) = qa(q) — r(q) (see [l(| E3 for details). In section III, we numerically analyze both 
r(q) and f(a) to detect the multifractality of the QKR-eigenfunctions. 



The Anderson model for a disordered system is described by a <i-dimensional disordered lattice, of size L, with 
a Hamiltonian H = X)n e ™ a n a ™ — S n #m ^ m ™( a n a ™ + a « a m) m the tight-binding approximation In the site 

representation, H turns out to be a sparse matrix of size N = L d with the diagonal elements as the site energies 
Hkk = £fc and the off-diagonals H mn = b mn given by the hopping conditions. For a Gaussian type on-site disorder 
(of variance u) and zero mean) and a nearest neighbor (n.n.) isotropic hopping with both random (Gaussian) and/ or 
non-random components, H can be modeled by an ensemble (later referred as the Anderson ensemble or AE) with 
following density 



C. Anderson Hamiltonian and the Complexity Parameter 



p{H, v, b) — Cexp 




(9) 



k 



(jfe,/)=n.n (k.l)^n.n. 
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with C as the normalization constant. As discussed in [19| . the above ensemble can be rewritten as 

P {H, h, b) = Ccxp[- 5^(l/2/ifci)(JT w - 6 fc; ) 2 ] (10) 



fc<; 



where = 0, ft-fefc = u>,hki — rjfu with /(fcZ) = 1 for {k,l} pairs connected by the hopping, f{kl) — > for all {fc,Z} 
pairs representing the disconnected sites. The single parameter A ae governing the spectral statistics (see eq.(19) of 
can then be given as 

with 

a w = ln|l- 7o w| + (2r/2)ln[|l- 270771^ + ^1] (12) 

with z as the number of the nearest-neighbors, on — — ln2 and 70 as an arbitrary constant. Further F(E) is the mean 
level density, C as the localization length and I 2 yp as the typical inverse participation ratio: oc (~ d . 

III. NUMERICAL ANALYSIS 

The objectives of our numerical analysis are two fold: (i) a search for the critical points of the quantum kicked 
rotor, and, (ii) a comparison of its fluctuation measures with those of a 3-dimensional Anderson ensemble. For this 
purpose, we analyze the following cases: 

(i) QKR1: quantum dynamics time-reversal but parity broken: k 2 = x 2 A, \ « 1.5,7 = 0,6*0 = w/2N,T = 
l,r = h = 40n/N which gives K ss 189/v^V- This case corresponds to the critical set (i) in section II. B and is 
analyzed for many sizes (N = 213 — * 1013) to verify the critical behavior. 

(ii) QKR2: quantum dynamics with both time-reversal and parity broken: k » 20000,7 = \ q N~ 3 ' 2 , \ = 
6,T = 1,6q = tt/2N,t — h — 40ir/N. This case belongs to the set (ii) in section II. B and its critical nature is also 
confirmed by analyzing many TV values. 



(iii) QKR3: : same as QKR1 but with x = 0.8 which gives K w 100/vAT. We consider this case to verify the 
analogy with a d = 3 Anderson ensemble. 

(iv) QKR4: quantum dynamics time-reversal but parity broken: k sw 4.5/?i,7 = 0,6*0 = (fiN~ 3 / 2 ,4> — 

0.847r 2 ,T = 1,t = h — 8ir/N (case (iii) in section II. B). This is also analogous to the Anderson system mentioned 
above in the QKR3 case, notwithstanding the crucial changes in k and 6*o for the two QKR cases. 

To explore critical behavior, we analyze large ensembles of the matrices U for both QKR1 and QKR2 for various 
matrix sizes N; the ensemble in each case is obtained by varying k in a small neighborhood while keeping N fixed. 
The chosen N range give K in the mixed regime for QKR1 (6.25 < K < 13) and in the chaotic regime for QKR2 
(1000 < K < 12000). Prior to the analysis, the quasi-energies (the eigenvalues of U) are unfolded by the local mean 
level density Z? -1 (= N/2ir, a constant due to repulsion and a unit-circle confinement of the quasi-energies). The 
figures 1,2 display the nearest-neighbor spacing distribution P(s) and the number variance S 2 (r), the measures of 
the short and long-range spectral correlations respectively, for the QKR1 and the QKR2. Note the curves in figure 1 
are intermediate to the Poisson and the COE limits; the size independence implies their survival in the infinite size 
limit too. This indicates the QKR1 as the critical point of transition from a localized phase to the delocalized phase. 
Similarly the curves in figure 2, intermediate to the COE and the CUE limits, suggest the QKR2 as the critical point 
of transition from the time-reversed phase to the time-irreversible phase (both phases in the chaotic regime). 

To reconfirm the critical nature of the QKR1 and the QKR2, we numerical analyze the moments of their local 
eigenfunctions intensity for various sizes. The results shown for r q in figures 3a, 4a indicate the multifractal nature 
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of the eigenfunctions; for small-<7 ranges, r q shows a behavior r q — (1 + c) q — d — c q 2 (or D q = 1 — cq) with d = 1 
and c = 0.06,0.075 for the QKR1 and the QKR2 respectively (see table 1 for the first few values of D q ). These 
results are reconfirmed by a numerical study of the /(a) spectrum displayed in figures 3b, 4b. For this purpose, we 
use the procedure based on the evaluation of moments, described in [l3| (using eq.(4) and eq.(10) of [T^D which has 
the advantage of full control over the finite-size corrections. The numerical results for a^,a\ and a\/2 for the QKR 
(see table 1) indicate a parabolic form of f(a) (also confirmed by the fits shown in figure 3b, 4b) and therefore a 
log-normal behavior of the local eigenfunction intensity u — |VK r )| 2 for large u-regions (see [ljl [HI EH f° r details). 
As shown in figures 3c, 4c, the tail-behavior of P u (u') can be well-approximated by the function 

f(u') = VW2^e bW e- a(u ' +c)2/2) . (13) 

with u' = [hut — (u)]/(ln 2 -u) (see figures 3c, 4c for numerical values of a,b,c). 

The bulk statistical behavior of the Hermitian matrices is known to be analogous to the unitary matrices (l?! |22| . 
This, along with the single parametric formulation of the statistics of the Hermitian matrix ensemble [l5j] , suggests the 
analogy of the QKR- Anderson ensemble statistics if their A parameters are equal (besides similar global constraints 
e.g. global symmetries) [l6j . Our next step is the comparison of the fluctuation measures of a time-reversal Anderson 
case with the QKR3, a time-reversal system with partially localized wavefunction in the momentum space (dynamical 
localization). For this purpose, we analyze a cubic (d — 3) Anderson lattice of linear size L (N — L d ) with a Gaussian 
site disorder (of variance u) — W 2 /12, W — 4.05 and mean zero), same for each site, an isotropic Gaussian hopping (of 
variance 77 = 1/12 and mean zero) between the nearest-neighbors with hard wall boundary conditions; these condition 
correspond to the critical point for a disorder driven metal-insulator transition [Hj|). A substitution of the above values 
(with t = 0) in eq. (Ti"2"]) gives a w — a, = 1.36. As shown by the numerical analysis in [n| [25|, F(E) ss 0.26e _£;2 / 5 and 
I^ vp ss 0.04 which on substitution in cq. pTj) gives A ae = 0.056 (with 70 = 2). 

For AE-QKR comparison, we analyze the ensembles of 2000 matrices with matrix size N = L = 512 for the 3d AE 
case and N = 513 for the QKR case. The energy dependence of A (see [19, 20]) forces us to confine our analysis to only 
10% of the levels near the band-center from each such matrix. The levels are unfolded by respective local mean-level 
density in each case (so as to compare the level-density fluctuations on a same density-scale) (l7j . Figure 5a shows 
the AE-QKR3-QKR4 comparison of P{s); the good agreement among the three curves verifies the A-dependence of 
the spectral correlations. This is reconfirmed by figure 5b showing a comparison of the number variance. Note that 
<j) ps 0.84-7T 2 for QKR4 is same as the theoretical analog given by the condition A kr ,t = A ae = 0.056 (see eas. (|6lll[) ). 
However x = 0.8 for QKR3 show a small deviation from the theoretically prediction (\ ~ 0.95). This may be due to 
eq.® being a poor approximation at the integrable-nonintegrable boundary K — 0. Note that the classical limits of 
QKR3 and QKR4 are different (K — for QKR3, K — 4.5 for QKR4), although their "semi-classical" K are same 
(K = 4.5 with N = 513 for QKR3, K = 4.5 for any N for QKR4). As clear from eq.©, K = marks the boundary 
between the integrable and mixed dynamics; K = 4.5 corresponds to the mixed nature of the dynamics (see figure 
5c). 

As discussed in (25|, the eigenfunction fluctuations of finite systems are influenced by two parameters, namely, 
system size iV as well as A measure . To compare A measure dependence of an eigenfunction measure, therefore, same 
system size should be taken for each system. Figure 6a shows the distribution P u (u') of the local eigenfunction 
intensity, u' = [hut — (lnu)]/(ln 2 i(), for the AE and the QKR3. The close proximity of the two curves suggests Ak r ,t 
as the parameter governing the local eigenfunctions intensity too; (we have verified the analogy also with QKR4). A 
comparison of the AE-QKR3 multifractality spectrum, shown in figure 6b, reconfirms their close similarity at least 
on the statistical grounds. 

The numerical confirmation of the statistical analogy of QKR4-AE-QKR3 systems supports our claim regarding 
single parametric (A)-dependence of the statistics besides global constraints. Note the latter are same for both 
QKR3 and QKR4 (i.e parity violated, time-reversal preserved and mixed dynamics) which results in their statistics 
intermediate to same universality classes, namely, Poisson and COE although the transition parameters are different 
in the two cases. 

The QKR-AE analogy can be utilized to connect them to other complex systems too. In [ll| [2f|, we studied the 
AE-connection with the PRBME (described by (H kl ) = 0, (H^) oc [1 + \k - l^/p 2 }' 1 [3) and the GBE ({H kl ) = 0, 
{H 2 t ) oc [1 + cN 2 ] -1 [l9|]). For the AE case considered above, the PRBME and the GBE analogs for the spectral 
statistics turned out to be p = 0.4 and c = 0.1; these systems are therefore the spectral statistical analogs of QKR3, 
QKR4 as well as of a N x N circular Brownian ensemble with w « 0.47rA -1 / 2 (see case(i) of section II. B). 
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IV. CONCLUSION 



To summarize, we studied the statistical analogy of two paradigmatic models of dynamical and disordered systems, 
namely, the quantum kicked rotor and the Anderson Hamiltonian, in the partially localized regimes. Our results 
indicate the existence of critical behavior in the classically mixed regime of the QKR with a smooth potential. This 
is qualitatively analogous to a disorder driven metal-insulator transition in the Anderson system; the quantitative 
analogy for their statistical behavior follows if their complexity parameters are equal. Our study also reveals the 
possibility of other transitions in the QKR e.g from a symmetry preserving phase to a symmetry fully violated phase. 
These transitions are analogous to similar symmetry breaking transitions in disordered metals e.g the Anderson 
Hamiltonian in the weak disorder limit in the presence of a slowly varying magnetic field. 

As with the Anderson transition, the QKR transitions are governed by the complexity parameter A too. However, 
contrary to the A-derivation for the Anderson case by an ensemble route, A for the QKR is derived by a semi-classical 
method. The semi-classical A-formulation is also numerically verified for the QKR-ensembles. This indicates an 
equivalence of the ensemble- averaging and the phase-space averaging for the statistical analysis. This further lends 
credence to the single parametric formulation of the statistical behavior of complex systems, irrespective of the origin 
of their complexity. However it needs to be examined for other dynamical systems. 

Research has indicated a multi-parametric dependence of the spectral-statistics at long energy scales of the dy- 
namical systems [26j, originating in the level-density oscillations due to short periodic orbits. However, these studies 
are not at variance with our work. This is because the "semi-classical" A-derivation in [20l | is based on the assumed 
equivalence of the traces of the operators with their phase-space averages. The assumption may not be valid on short 
time-scales of the dynamics. One should also understand the exact role of the ensemble-averaging for the statistical 
analysis of dynamical systems. 
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VI. FIGURE CAPTION 

Fig. 1 [Color online]. Spectral Measures of QKR1: 

(a) the nearest-neighbor spacing distribution P(S) of the eigenvalues for various sizes N, with inset showing behavior 
on the linear scale, (b) Number variance S 2 (y) for various sizes. 

The convergence of the curves for different sizes indicates scale-invariance of the statistics. The behavior is critical, 
being different from the two end-points, namely, Poisson and CUE statistics even in infinite size limit. 



Fig. 2 [Color online]. Spectral Measures of QKR2: (a) Distribution P(S) of the nearest-neighbor eigenvalue spacings 
S for various sizes, with inset showing behavior on the linear scale, (b) Number variance £ 2 (r) for various sizes. 

Again the statistics being intermediate between COE and CUE, and, convergence of the curves for different sizes 
indicates its critical behavior. 



Fig. 3 [Color online]. Multifractality of QKR1 : 

(a) Fractal Dimension r q along with the fit y(q) = (1 + c)q — 1 — cq 2 with c = 0.06 (good only for q < 3). A fit for the 
large q regime suggest following behavior: r q = q — 1 + 0.02q 2 . (b) Multifractal spectrum f(a) for various sizes along 

with the parabolic fit f(a) — d~ fe^j^ wrfcn a o = 1-09 and d = 1. (c) Distribution P u {u') with u' = [hm— (u)]/(ln 2 w) 
of the local intensity of an eigenfunction for QKR1. The solid line represent the function f(u') given by eq. (fT3")) with 
a = 5.2, b = 1.2 and c = 0.78 (corresponding to an approximate log-normal behavior of P u (u)), a good approximation 
in tail-region as expected. The inset shows the behavior on linear scale. 



Fig. 4 [Color on line]. Multifractality of QKR2: (a) Fractal Dimension r q along with the fit y(q) = (1 + c)q — 1 — cq 2 
with c = 0.075 (good only for q < 4). (b) Mulifractal spectrum f(a) for various sizes along with a parabolic fit, of the 
same form as in figure 3.b with ao = 1.045 and d = 1. (c) Distribution P u (v!) of the local intensity of an eigenfunction 
for QKR2, with u' same as in figure 3c. The solid line represent the function f(u') with a = 5.3, b — 0.95 and c = 0.73 
(corresponding to a log-normal behavior of P u (u)) which fits well in the tail region of P u (v!) . The inset shows the 
behavior on the linear scale. 



Fig. 5 [Color online]. Comparison of spectral statistics of the Anderson ensemble with QKR3 and QKR4: (a) P(S), 
with inset showing the linear behavior, (b) X 2 (r). Note QKR3 and QKR4 turns out to be close to the one suggested 
by the relation A a = Ak r ,t (giving \ — 0-95 for QKR3 and 4> — 0.84-7T 2 for QKR4). (c) Phase-space behavior of the 
classical kicked rotor at K — 4.5 (see eq.([2])) 



Fig. 6 [Color online]. Comparison of the multifractality of the eigenfunctions of QKR3 with the Anderson ensemble: 
(a) local intensity distribution P u (u') of an eigenfunction with vl same as in figure 3,4. Also shown is the function 
f(v!) with a = 15.9, b = 7.7, c = 0.99, a good fit in the tail region of P u {u'). The inset shows the behavior on linear 
scale, (b) multifractal spectrum f(a): Note the QKR3 analog of the Anderson ensemble here is same as the one in 
figure 5. Again the fit has the same parabolic form as in figures 3,4, with ao — 1.1 and d = 1. 
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TABLE I: Multifractality Analysis of QKR: a values here are obtained by an L — > oo extrapolation of The < T)% < 1- 

behavior indicates a multifractal nature of the three QKR cases. 
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